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If left-right gauge theory, SU{2)l x SU(2)r x U{1)b-l x SU{3)c{gL = gR.){= G2213), occurs as 
an intermediate symmetry in a grand unified theory then, apart from other advantages, it is possible 
to obtain the seesaw scale necessary to understand small neutrino masses with Majorana coupling 
of order unity. Barring threshold or non-renormalizable gravitational effects at the GUT scale, or 
the assumed presence of additional light scalar particles of unprescribed origin, all other attempts 
to achieve manifest one-loop gauge coupling unification in SUSY 5*0(10) with such intermediate 
symmetry have not been successful so far. Attributing this failure to lack of flavor symmetry in the 
grand unified theory, we show how the spontaneous symmetry breaking of 50(10) x S4 leads to such 
left-right intermediate breaking scale extending over the range Mfl ~ 5 X 10^ GeV to 10^^ GeV. 
All the charged fermion masses are fitted at the intermediate see-saw scale, Mn — Mr ~ 4 x 10^^ 
GeV which is obtained with Majorana coupling /o — 1. Using type I seesaw and a constrained 
parametrisation in which CP-violation originates only from the quark sector, in addition to other 
predictions made in the neutrino sector, the reactor mixing angle is found to be S13 ~ 3° — 5° 
which is in the range accessible to ongoing experiments. The leptonic Dirac phase turns out to be 
S ~ 2.9 — 3.1 radians with the predicted values of Jarlskog invariant Jcp — 2.95 x 10"''' — 10""^. 

PACS numbers: 14.60.Pq, ll.30.Hv, 12. 10. Dm 



I. INTRODUCTION 

SO{10) grand unified theory [l| has a number of attractive features which have resulted in recent surge of inves- 
tigations including applications to fermion masses and mixings 2]. It unifies all fermions of one generation plus the 
right-handed neutrino into one spinorial representation 16. With D-parity as an element of gauge transformation [sl, 
it naturally restores left-right and CP symmetries at the GUT scale and thus, it can provide a spontaneous origin of P 
and CP-violations Q ■ It embodies quark-lepton unification with high predictive power in the fermionic sector [^, 0] . 
Through its Higgs representations 10 and 126 or 16, it has the potentialities for intermediate SU{2)ii x U{1)b-l 
breaking, generation of small Majorana neutrino masses through type-I, type-II, and type III see-saw mechanisms 
0i B @| J explanation of large neutrino mixings through type-II see-saw dominance, accounting for dominant contribu- 
tions to charged fermion masses through the 10-representation and providing the desired corrections to them through 
the weak-doublets in 126 and 120 [ifl, [H [H, [3, Q [H, [B . With natural R-parity conservation, in addition to 



ensuring proton stability [17i |. it predicts the lightest supersymmetric particle as a stable dark matter candidate. 

All fermion masses and mixings, including very small masses and large mixings in the neutrino sector, have been 
shown to fit reasonably well if the right-handed neutrino mass scale is in the range 10^'^ GeV to 10^^ GeV [l^, ■ 
Further, thermal leptogenesis explaining origin of matter through baryogenesis can be implemented if the right handed 
neutrino masses are in intermediate range [isl . [l9| . Therefore, it would be interesting to obtain the right-handed mass 
scale near the intermediate breaking of SU{2)l x SU{2)ji x S't/(4)c(= G224) C 5*0(10) or any of its subgroups like 
G2213- However, the mass spectra analysis in the minimal SUSY 5*0(10) model with the Higgs representations 
210 © 126 ® 126 © 10 rules out any possibility of intermediategauge symmetry by predicting additional light scalars 
which are found to disrupt gauge coupling unification 21, 2^, 2^. It has been also noted that, even if these additional 



light scalars are made naturally superheavy using a non-minimal Higgs representation, leaving only the minimal 
light Higgs spectrum necessary to implement spontaneous symmetry breaking in the presence of supersymmetry and 
R-parity conservation, there exists no intermediate scale through manifest one-loop unification of gauge couplings. 

Recently attempts have been made to obtain desired see-saw scales and improved fits to the fermion masses using 
different types of mechanisms or by invoking GUTs in higher dimensions 15|, |24|. In addition, extensions of gauge 
symmetries by non-abelian fiavor groups like 53,54, and have resulted in interesting consequences 0,11^ 2^, 27, 
28 . 2^ including high scale unification of quark and lepton mixings [ij- 

Since SUSY GUTs with a left-right intermediate gauge symmetry has many attractive features over super-grand 
desert models, in this paper we discuss a novel procedure of constructing such a model through an extension of the 
left-right gauge symmtry to G2213 x 54 and the corresponding extension of the GUT symmetry to 50(10) x 54 to 
encompass supersymmetric flavor unification, grand unification, and R-parity conservation (30| . We find that when 
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such a flavor symmetry is included, manifestly successful one-loop gauge coupling unification occurs with the three 
gauge couplings of the left-right gauge theory attaining convergence to the unification coupling at the GUT scale. 
The intermediate scale is predicted over a wide range with « 5 x 10^ GeV to 10^^ GeV. The desired Higgs 
scalars necessary for the gauge coupling unification are found to be consistent with the mass spectra analysis in the 
S'O(IO) X Si model. Although the contribution of three fermion generations cancel out from one-loop unification, 
flavor symmetry requires enlarged Higgs spectrum which modifles the beta function coefficients and the evolution of 
the gauge couplings leading to successful unification in the presence of x 'S'4 intermediate symmetry. In the 

second part of the paper, we fit the renormalization group (RG) extrapolated data on fermion masses, mixings, and 
phases at the intermediate see-saw scale Mm ~ Mr ~ lO^'^ GeV and obtain successful predictions in the neutrino 
sector. 

This paper is organized in the following manner. In Sec. 2 we briefly review problems associated with realization 
of an intermediate scale in SUSY 50(10). In Sec. Ill we discuss how the desired intermediate scale is achieved via 
left-right gauge and 84^ flavor symmetries. Symmetry breaking of 5*0(10) x ^4 is discussed in Sec. IV. In Sec.V we 
show how the desired light particle spectrum is obtained from the 50(10) x 54 theory. Fits to the fermion masses 
and mixings and model predictions in the neutrino sector are carried out in Sec. VI. A brief summary of investigations 
made and conclusions obtained are stated in Sec. VII. 



II. DIFFICULTIES IN R-PARITY CONSERVING LEFT-RIGHT INTERMEDIATE GAUGE SYMMETRY 

In this section we discuss briefly the problems associated with obtaining a SU{2)l x SU{2)ii x U{1)b-l x 
SU{3)cigL — 5i?)(= G2213) intermediate gauge symmetry in R-parity conserving supersymmetric SO(IO) grand 
unified theory through one-loop unification of gauge couplings. 

For this purpose, we consider two-step breaking of 50(10) to the minimal supersymmetric standard model(MSSM) 
through G2213 intermediate gauge symmetry in the so called minimal grand unified theory. 



21 n 

50(10) — G2213 

Mu 

G213 U{1U X SU{3)c. (1) 

Mft Mz 

The G224 sub-multiplet (1,1,15) in $(210) contains a G2213 singlet which is even under D-parity Q. When this 
component acquires VEV, 50(10) — > G2213 with unbroken left-right discrete symmetry. 

Unlike the D-parity breaking case where the intermediate left-right gauge group has four different coupling constants, 
in the present case G2213 has only three gauge couplings, (72L = 52fi,<73C, and gBL for [i > Mu. In the second step, 
the right-haded triplet component in 126 acquires VEV to break G2213 G213 while generating heavy right-handed 
Majorana neutrino mass. In the process of spontaneous electro-weak symmetry breaking driven by weak bi-doublet 
in 10, small left-handed neutrino masses are generated through Type I and Type II see-saw mechanisms. 



To discuss gauge coupling unification we use the following three RGEs from Mz to Mjj |31l. |32|. |3 



1 1 ay Mr 1 . , , . Mjj 

—TW-\ = — + T^1^7/ + T7r 3°2L + 2a'Bi In-/, 
aY{Mz) OiG 2t: Mz IOtt Mr 

1 = J_ I In + ^ In — 
a2L{Mz) ^ oiG 2tt ^ Mz 2tt ^ Mr' 

l_ ^ ^^a^^^MR_^a^^^Mu^ 
a3c{Mz) ac 27r Mz 2-n Mr 

where ac is the GUT fine-structure constant and the beta function coefficients and are determined by the 
particle spectrum in the ranges from My to M^, and from Mr to Af[/, respectively. Adopting the standard procedure 



we obtain the following two equations [3 



27r / % S^n'ew(.Mz)\ ,,M^ , „,„M. 



U ^ -7^(1- ""^-^7 ^"'^M = ^In^ + Bln^, (3) 



a(JMz) \ 3 y Mz M, 

27r / 8a(Mz) \ _ A% , „,,„Mk 



Ls = -^^[l- = A'ln^ + S'ln^. (4) 

a{Mz) V iazciMz)) Mz Mz ^ ' 
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where 



^ = ^(.<^'bl- a^L), B =^{aY - a2L) -^{a'sL- a2L), 



2 8 5 8 2 o 

A' = 2a'2i + -a'sL - -^a'^c^ B' = -ay + - -age - (2^^ + -a^^ - -4^). (5) 



From eq.([3]) and eq.(j4|), the analytic expressions for Mu and Mr immediately follow, 



Using PDG values, a{Mz) = 127.9, S'm26'vi/(A^z) = 0.2312 , and a3c{Mz) = 0.1187 [3J, we obtain 



662.736, Lg = 308.305. (8) 

Ignoring the contributions from additional lighter scalar multiplets emerging from mass spectra analysis which have 
been discussed later, the minimal Higgs content necessary to break SUSY S'O(IO) through G2213 to the low energy 
group and the associated beta function coefficients are, 
/i = Mz - Mr: 



H"(2,l,l)®i/'^(2,-l,l) C G213, 
ay = 33/5, a2L = 1, asc = -3, 



(9) 
(10) 



/i = Mr - Ml, 



H^{2, 2, 0, 1), Al(3, 1, -2, 1) ® Afl(l, 3, -2, 1) ® Al(3, 1, 2, 1) ® Ar{3, 1, 2, 1) 



(11) 



under Go 



with 



^BL — 24, — a2R — 5 

Thus, the one-loop coefficients give 



' - aac = -3. (12) 



A = 38/3, B = -10/3, A' = 34, B' = -14, AB' - A'B = -64, (13) 
leading to the solutions [s^, 

Mr ~ 10^^ GeV, Mu = 2x 10^^ GeK (14) 

The mass spectra analysis in the minimal SUSY 5*0(10) with 210® 126® 126® 10 predicts that there are additional 
scalar components of 210 having intermediate scale mass with the following G2213 quantum numbers [23| . 

(3, 1, -2/3, 3), (3, 1, 2/3,3), (1, 3, -2/3, 3), (1, 3, 2/3, 3) (15) 

Not onl y th ese states prevent any value of intermediate scale below Mu, as explicitly noted in Ref. [23], but, as noted 
in Ref. [35], their presence at scales substantially lower than 2 x 10^® GeV also spoils perturbative renormalization of 
gauge couplings. 
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It has been pointed out that G2213 intermediate scale can still be obtained in non-minimal SOCIO) by threshold 
effects or by the presence of non-renormalizable dim. 5 operators in the SO{10) Lagrangian 35|, |3g|. These might 



arise if, in addition to 210, the theory contains a Higgs representation 54. The implementation of the gauge coupling 
unification has been found to be possible by threshold corrections or by gravitational corrections in non-renormalizable 
SUSY SO(IO) only if the additional light scalars given in eq. are made superheavy which could be realized due 
to the added presence of 54 [35]. 

However, the more attractive and popular unification scheme being through manifestly one-loop evolution of gauge 
couplings in a renormalizable grand unified theory where coupling constants from lower scale evolve to converge to the 
unification coupling at the GUT scale, in the next two sections we show how such a unification is achieved when the 
flavor symmetry ^4 is combined with SO{10) as well as the R-parity and parity conserving supersymmetric left-right 
gauge theory G2213 at the intermediate scale. 

III. INTERMEDIATE LEFT-RIGHT GAUGE SYMMETRY WITH S4 FLAVOR SYMMETRY 

In this section we show that in the presence of flavor symmetry and left-right gauge symmetry, the extrapolation of 
standard model gauge couplings through G2213 x S4 intermediate symmetry naturally leads to successful unification 
of gauge couplings at 2 x 10^^ GeV. In the next section we show how the minimal particle content necessary for this 
intermediate symmetry follows from mass spectrum analysis of supersymmetric S'O(IO) x S'4. 

It is well known that the minimal particle content of the standard model (SM) alone in non-supersymmetric theory 
does not allow its gauge couplings to unify at any higher scale. However, when the SM emerges from nonSUSY 
left-right symmetric G2213 at the intermediate scale, profound unification of gauge coplings occurs at the GUT scale 
[371] . Since the non-SUSY theories are well known for their generic gauge hierarchy problem, when supersymmetry 
is combined with the SM to solve this problem, the enlarged particle spectrum naturally infused into the MSSM, 
achieves gauge coupling unification at the SUSY GUT scale without any intermediate scale. On the other hand, 
if an intermediate symmetry with extended gauge group like G2213 is introduced, manifest one-loop gauge coupling 
unification is spoiled in SUSY GUTs. 

Earlier, without ascribing any connection with flavor symmetry, several authors have noted that when the par- 
ticle spectrum at lower scales is further extended beyond the minimal spectrum, manifest one-loop gaug e coupling 
unification occurs in the presence of G2213 intermediate symmetry with or without parity or R-parity 



These observations lead us to suggest that the present failure to achieve manifest one-loop gauge coupling unification 
with intermediate left-right gauge symmetry may be hinting at its extension to include a family symmetry with 
corresponding extesion in the particle spectrum. We find that this new symmetry to be appended to G2213 and 
50(10) could be the weU known flavor symmetry S4 0, [13, 111]. 

This conclusion can be inferred by also looking into the structure of the RGEs in eq.Q and eq.(|4]). It is clear that 
if the coefficients B and B' are negligible compared to A and A', then the solutions for the mass scales would be 
insensitive to the values of M^. Then values of Mn substantially lower than the SUSY GUT scale could be tolerated 
by the RG constraints. We find that this possibility can be realised within 5*0(10) x 54. 

In the enlarged particle spectrum, along with the two sets of Higgs triplets of the minimal scenario given in ea. (fT2)) . 
the presence of 54 symmetry needs six bi-doublets instead of only one [28l |. In additon, other scalar multiplets 
belonging to 1 , 2 , or 3 of 54 having nontrivial transformation property under G2213 are also found to be essential. 
These latter Higgs particles turn out to be a triplet 3 of 54 and each member of the trplet transforms as a color octet. 
The fermions of three generations are taken as 3' of 54. 

Thus, keeping the MSSM particle spectrum from Mz to Mn unaltered, the enlarged Higgs spectrum at the 
intermediate scale consistent with G2213 x 54 symmetry is: 

11 = Mr- Mu: 



Al(3, 1, -2, 1) ® Afl(l, 3, -2, 1) ® Ai(3, 1, 2, 1) ® Afl(3, 1, 2, 1), 

6(2,2,0,1), 3(1,1,0,8), (16) 



where 6 = 3 + 2 + 1, and 3, 2 and 1 are triplet, doublet, and singlet , respectively, under 54. These Higgs scalars 
modify the beta-function coefficients to 
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= 24, a'2L = a'2R = 10, 4^ = 6 



(17) 



Noting that the value of a'g^r = 24 which is the same as in eg . (fT2|) . determines the perturbative constraint on the 
lowest allowed value of Mn [35j , we have with the above Higgs spectrum. 



28 

A^—, A' ^20, B= B' = 0. 



(18) 



It is interesting to note that the particular combination of G2213 x 6*4 given in ea. (|16p leads to exactly vanishing 
values of B and B' . The fact that now B and B' both vanish with such Higgs content ensures the possibility of an 
intermediate scale over wide range of values. But the perturbative lower bound on Mji being determined by a'gj^ = 24 
due to the appearance of a Landau pole, now the popular one loo p u nification of gauge couplings is expected to 
materialise for values of the intermediate scale satisfying this bound (35| . 

As eq.® and eq.(I7]) are no longer valid with B — B' = , we solve for the mass scales numerically using eqs.(2), 
(10), and (17). We find all values of the left-right symmetry breaking scale are permitted over a wide range. 



5 X 10^ GeV < A/r, < 10^^ GeV 



(19) 



but having almost the same value of unification scale Mjj = 2 x 10^® GeV for all solutions. Two examples of such 
solutions for — 10^^ GeV and M/? = 5 x 10^ GeV are shown in Fig. [1] and Fig. [21 respectively, where manifest 
one-loop unification with the three gauge couplings of G2213 converging at the GUT scale, Mjj = 2 x 10^^ GeV, is 
evident. 



U(l)v 


Mjj ^10'-^ GeV 








iU(3)c 













FIG. 1: Evolution of gauge couplings showing variation of inverse fine-structure constants as a function of the mass 

scale fi with G2213 x 54 intermediate gauge symmetry breaking at Mr = 10^^ GeV in SUSY SO(IO) x 5*4 model. The top solid 
line represents ay^(/i) for fi = Mz — Mr and ag^_j^{^) for n = Mr — Mu- The middle and the bottom lines represent a^^(/x) 
and a^Q{fi), respectively, throughout the range of /i. 



For Mfj ~ 10""^^ GeV , the GUT fine structure constant has the value aa — 1/24.5 which increases as approaches 
lower values, reaching ac — 1/12 and ac — 1/2 at Mr — 10^^ GeV and Mr = 5 x 10^ GeV, respectively. This 
phenomenon is due to the appearance of Landau pole near Mu ~ 2 x 10^^ GeV in the gauge couphng of U{1)b-l 
when Mr < 10^ GeV. Now that all the three gauge couplings of G2213 are unified at the GUT scale, ac would 
approach 00 for the same value of Mr < 10^ GeV. The basic reason is that a'g_j^ = 24 has remained the same as the 
minimal model inspite of new contributions from G2213 x ^4 multiplets. 

The SU{2)r gauge couphng at ^ ~ Mr is nearly gi^ — qr = g :^ 0.7. With the SU{2)r gauge boson mass 
scale Mr ~ gvR, the vacuum expectation value of the Higgs triplet field is also allowed in the similar range with 
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FIG. 2: Same as Fig. [T]but for Mr = 5 x 10^ GeV. 



TABLE I: Particle content of the model and their transformation properties under S4 x 50(10) 



Fermions 




Higgs Bosons 












*i,(i = l,2,3) 


S 




Al_2,3 


So e So 


Ho 


Ul,2 


H3,4,5 


3' X 16 


1 X 54 


1 X 210 


3 X 45 


1 X 126 126 


1 X 10 


2 X 10 


3 X 10 



< Ali >— ~ 5 X 10^ GeV to 10^^ GeV. Ignoring the constraint from the neutrino osciUation data which wiU be 
discussed in Sec. 6 , the right-handed neutrino mass Mn ~ faVR for Majorana coupUng /o ~ 1 is then aUowed to vary 
over similar range without the necessity of any tuning of /q. 

In the next section we show how the minimal particle content needed for gauge coupling and flavor unification 
through 6*2213 x *S'4 can be easily embedded in 5*0(10) x 6*4. We further show how the F-term flatness condition leads 
to the desired minimal particle spectrum of the model below the GUT scale while keeping all other components of 
the 5*0(10) X Si representations superheavy. 



IV. SYMMETRY BREAKING OF 5*0(10) x 5*4 
We consider 5*4 flavor symmetry for three fermion generations and supersymmetric grand unification of three forces 



of nature through 5*0(10) x 54 [27|, |28|. Instead of starting from flavor symmetric standard model without SUSY 
discussed in Ref. ^2811, we consider flavor symmetry starting at the intermediate scale through G2213 x ^4. Such a 
model leading to MSSM at lower scales originates from intermediate breaking of 50(10) x 54, 

50(10) X 54 G2213 X 54 Gsis^t f/(l)e.. X 5t/(3)c. 

To achieve successful gauge coupling unification with the desired intermediate scale, we need the minimal particle 
content same as in the MSSM from Mz to Mji. In the presence of G2213 x 54, we need six bi-doublets each of which is 
embeddable in _ff (10)o, -ff (10)i.2, i?(10)3.4^5 of 50(10) and these transform as singlet, doublet, and triplet, respectively, 
under 5*4. The left and the right-handed Higgs triplets needed to maintain supersymmetry and implement spontaneous 
breaking of G2213 at scale Mr are contained in Eo(126) and So(126) of 50(10). We note that, in addition to the 
minimal number of doublets and triplets, the gauge coupling unification also needs three 5C/(3)c octets transforming 
as (1,1,0,8) under G2213. Each of them is contained in the G224 multiplet (1,1,15) which, in turn, is contained in 
the Ai{Ab) (i — 1,2,3) of 50(10) treated as a triplet 3 of 54. All other Higgs particles are to become superheavy 
with masses near the GUT scale for successful gauge coupling unification at one-loop level and this can be achieved 
in the presence of 5(54) of 50(10) [ij]. The Fermion and Higgs representations are given in Table HI 

In order to realize such a spectrum by actual potential minimization in the presence of supersymmetry, we break 
50(10) by giving GUT-scale vacuum expectation values to the two D-parity conserving singlets of $(210) and 
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S'(54). We follow the mass spectrum analysis technique for supersymmetric S'O(IO) grand unification [2l|, 0, [2^ . 
We assign nearly equal vacuum expectation values to both the singlets with {S) ~ (<i>) such that, effectively, the GUT 
symmetry breaking to G2213 appears as one step process. 

Although for the sake of gauge coupling unification alone, it is possible to treat the effective gauge symmetry below 
M/{ as MSSMxS'4 with only one weak bi-doublet having mass near the electro- weak scale, for accommodating all 
fermion masses and mixings in this model as dicussed in the next section, it is necessary to break the flavor symmetry 
at the intermediate scale. 

We assume that 02213 gauge symmetry is broken at the scale Mn due to vacuum expectation value of the standard 
model singlet contained in the right handed triplet Ai^(l,3, 2,1) in 126. After the G2213 x S4 breaking, only two 
MSSM Higgs doublets are taken to remain light. 

The superpotential near the GUT scale can be written as, 

Wh = ^TO<E.$' + ^msS^ + V 4' + msSoSo + \mHoHl 
III — 'i I 

+S{X5S^ + As V A' + At*' + Xs^l + AgS^ + Aio-ffo' + ^iiH^ + X^H^). (20) 
^ — 

Using vacuum expectation values < S >= vs,< $ >= vq,,< cr, < >— a, the vanishing F-terms yield 




Iqv^ liaa 2ljv<j,vs 



3\/2 10\/2 Vis 

2V5 VTs 



msvs 



ms + 



IOV2. 



0, 

0, 
0. 



(21) 
(22) 
(23) 



Due to the vanishing D-term, a — a = vr and the corresponding F-terms for a ot a yield the same equation as 
eq. (f23|) . In the desired hierarchial case, both a and a are much smaller compared to < 5* >, < $ > leading to the 
relation between the GUT-scale VEVs and m$. 



3V2 VTs 

Using from ea. p4)) in eg. ([22)1 gives a quadratic equation for vs, 



0. 



(24) 



pvg + qvs - r = 0, 



where 



2V5 



ms 



24/f 
^ 5Vl5/g' 
241^ 



18/t 



15/^ 



(25) 



(26) 



In the next section we discuss the emergence of mass spectra necessary to keep only the desired minimal number 
of Higgs particles light while making others superheavy. 
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V. LIGHT AND HEAVY PARTICLE STATES FROM MASS SPECTRA 



In this section we discuss the emerging mass spectra from the spontaneously broken flavor symmetric GUT while 
making provisions for would be Goldstone bosons and the light scalars necessary for gauge coupling unification. In 
contrast to the minimal model without flavor symmetry where unwanted light scalar degrees of freedom are found to 
spoil gauge coupling unification (23| . in the present case, due to the presence of the scalar multiplet 54 in the Higgs 
superpotential, it is possible to lift those masses to the GUT scale. 



A. Goldstone Bosons 



In the process of spontaneous symmetry breaking of 50(10) ^ G2213 through < S > and < $ >, 30 gauge bosons 
would aquire GUT-scale mass by absorbing the corresponding mass-less scalars. Under G2213 these superheavy gauge 
bosons have the quantum numbers (2, 2, 2/3, 3) + (2, 2, -2/3, 3) + (1, 1, 2/3, 3) + (1, 1, -2/3, 3). Whereas the first two 
sets of states are contained in both (2, 2, 6) C G224 C 54 and (2, 2, 10) C G224 C 210 of 5*0(10), the next two sets of 
states are contained in (1, 1, 15) C G224 C 210 or 45 of 50(10). Using the superpotential in eq. (pn|) and the vacuum 
expectation values, we show how the desired Goldstone bosons are obtained. 



A.l. (1, 1, 3, 2/3) + (c.c) as Goldstone Bosons 



Noting that a — a = vr << vs ^ ^ Mjj, it turns out that these unmixed states in the leading approximation 
have masses. 

Using eq. ([24]) it is immediately recognised that these are naturally the light pseudo Golstone bosons to be absorbed 
by the mass-less vector bosons to make them superheavy. It can be easily checked that other unmixed states having 
the same quantum nubers in Ai{i — 1,2,3) have degenerate superheavy masses. 



which are naturally near the GUT scale. 



A. 2. (2, 2, 3, 1/3)+ (c. c) as Goldstone Bosons 

TT • 1 ■ If A \(2.2.1/3,3) /. n(2,2,1/3,3) n(2.2.1/3,3) £,(2,2,1/3,3) -^(2.2.1/3,3) ^ (2.2.1/3,3)i TZ^i , 

Usmg the basis [(Ai)^^^^^^! ■ ' , iA2)l^[2]e'j . (A3)(2.2,6) ' ^^{2,2,6} ' *(2,2,6) ' ' *(2,2:io) ] ^22] where the 
superscripts(subscripts) refer to gauge quantum numbers of the Higgs multiplets under 02213(0224) , there are three 
unmixed pairs of states in Ai and A2 and A^ with degenerate superheavy masses. 



niA 



2V15 



Vs. 



The fourth unmixed pair is that of ^'(2'2'g)''^''^^ with superheavy mass. 



7I7 



The remaining two pairs of states, 



Ml 



ms 



TO<J -I r= 

4^/15 
(c.c.) and $ 



VS. 



(2,2,1/3,3) 
(2,2,10) 



(29) 



(30) 



2V3 



2V5 



™* + 3V2 



(c.c), mix through the mass matrix, 



4%/5 



(31) 



It is clear that one linear combination of these two pairs can be made mass-less by tuning the parameter such that 
it supplies the remaining Goldstone modes. The other orthogonal combination acquires mass near the GUT scale. 
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B. Light Scalrs from Mass Spectra for Gauge Coupling Unification 



In Sec. 3, the unification of gauge couplings with G2213 x S4 intermediate gauge symmetry has been shown to 
require the usual left- and the right-handed triplets that are contained in 126 126, six bi-doublets contained in 
six 10-plets, and a set of three color octets transforming as (1,1,0,8) under G'22i3- These octets are contained in 
the 6*224 submultiplet (1,1,15) of ^^45) or $(210) of 50(10). In addition, the G224 submultiplet (1,1,20') in 
5(54) also contains the octet component. But we will find it covenient to obtain these three octets from a triplet of 
45i(? = 1,2,3) C 50(10). 

At first it is to be noted that the triplets Ai(3, 1, —2, 1), A/j(l, 3, —2, 1) and their conjugates contained 126 ® 126 
acquire degenerate masses. 

The condition Mji << Mfj can be ensured by tuning li. At first the six bi-doublets from the six 10-plets are 
treated to have masses near in the usual fashion by some doublet triplet splitting mechanism or by tuning the 
parameters lio, In, and I12 while the weak bi-doublets in 126 © 126 and 210 are kept heavy at the GUT scale. In the 
next section we show how five linear combinations of these bi-doublets can be treated to have masses at the M^^-scale 
while keeping the mass of the remaining linear combination at the electro-weak scale, thus supplying the pair of two 
MSSM doublets (iJ",77''). 

Choosing the basis iAi)i'i'^f,Sl'l'2Q'T^i'i'i5j ^^'^ that there are three unmixed states in Ai with masses. 



ITIA 



^/2l4V<!> 21qVs 
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(33) 



Clearly the advatage of Ai being the members of 3 C S4 is that the tuning of the single parameter I4 makes all the 
three octets light having masses near the AIji scale which is essential for gauge coupling unification. It is found that 
the other two states mix through the mass matrix, 



Mo = 



ms 



2\/3hvs 



3V2 Vl5 



(34) 



The eigen values emerging from eq. p4p are at the GUT scale and we do not adopt any further fine-tuning. 

We have verified that all the components from G224 multiplets,(l, 3, 15) © (3, 1, 15) C 210 acquire masses near the 
GUT scale due to the presence of 54 in the model and there are no other lighter states which are likely to disrupt 
gauge coupling unification. 

In summary the theory has enough parameter space for the successful implemetation of the model with G2213 x 54 
intermediate symmetry and gauge coupling unification at the GUT scale via 50(10) x 54. 



VI. FERMION MASSES AND MIXINGS 



In this section we address the question of fermion masses and mixings and make predictions in the neutrino sector. 
For this purpopose we assume type-I seesaw dominance and also include a pair of 126j_.2 © 126i^2 = ^1,2 © ^1,2 as 
members of an 54 doublet but with all their components having GUT scale masses [32] . Even though all the weak 

bi-doublets A; (i — 0, 1, 2) contained in the G224 submultiplets (2, 2, 15) of Si, {i — 0, 1, 2) have GUT scale masses, it 
is necessary to clarify how their induced VEVs contribute to the Dirac masses of all fermions in a manner analogous 
to minimal SUSY 50(10) 0. 



A. Light Weak Bi-doublets and Vacuum Expectation Values 

In this subsection we clarify how by keeping only the desired minimal number of particles below the GUT scale, 
VEV of different weak bi-doublets are made available to generate fermion masses withot disrupting gauge coupling 
unification. 

The added 54— doublet fields Sd © will have a new contributions to the superpotential which include, 
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W'h = ms'S^Sc + +\[<^>^l^D + {\2^l + X's^DHd^ (35) 

where ellipses denote couplings to other Higgs representations and subscript D stands to indicate a £'4— doublet. Since 
TOs' — Mij and we need no fine-tuning of parameters X'l 2 3 to maintain gauge hierarchy, it immediately follows that 
all the components of the S'4— doublet pair (including the weak bi-doublets and triplets) acquire masses at the GUT 
scale and they do not upset gauge coupling unification. Denoting the full superpotential as W — Wh + W^, the 
F— term due to # now contributes the following terms to the scalar potential, 

V — (A2S0 + A3Eo)ffo5^0^0 

+ {X'2^l + X'3^)Hd^o^o- (36) 

This has the implication that whenever the RH-triplcts in Sq © Sq and the weak doublet in S4-singlet Hq C IOq 
aquire VEVs vr and = yo < Hq > , respectively, where the latter is approximately of the order of the weak scale, 
the weak doublets in Sq gets an induced VEV. On the otherhand, in addition to vr, the VEVs of the S4— doublet 
components Hi 2 C 10i^2 generate induced VEVs in the weak doublets in Si_2- The order of magnitudes of all the 
three induced VEVs can now be expressed as, 

(Af) = ^^,(*-0,1,2). (37) 

where, as defined subsequently in this section, stands for the product of i*'' VEV and the respective Yukawa 
coupling m . 

One major difference from the minimal SUSY 50(10) is that in addition to the VEVs of up- and down type doublets 
of Ho, the VEVs of nontrivial S4— representations also enter into the RHS of eq. ([37]) . 

The next point that needs explanation is how only the six bi-doublets lighter than Mu acquire nearly electroweak- 
scale VEVs while five of them have masses near Mr scale and only the remaining bi-double has mass near the 
electro-weak scale to supply the up-type and the down-type MSSM doublets {H^,H'^). 

In order to achieve this objective we introduce two S'O(10)-singlet scalar fields, 775 and 77^ which transform as 
doublet and triplet, respectively, under 5*4. These will make additional contribution to the superpotential, 



W'j^^XsVsHQHD + X'srfsHoHT + .... (38) 

We assign order Mr scale VEVs to rjg and T]g to break the S'4 symmetry at the intermediate scale and generate 
Hq — H]j and Hq — Ht mixings. Then using bi-unitary transformation on the six doublet fielda to diagonalize the 
bi-doublet mass matrix at the intermediate scale, we treat only one linear combinations of the weak bi-doublets 
to have mass at the electro-weak scale while the remaining five linear combinations are treated to acquire mass at 
the intermediate scale. This would involve only one fine-tuning of the new parameters of the superpotential. Once 
the lightest linear combination of all the six bi-doublets is constructed in this manner to supply the MSSM Higgs 
doublets (iJ",7f'^), the electroweak VEVs of the latter imply VEVs of approximately the same order for all the bi- 
doublet components in i?j(i = 0, 1, 5). Then the induced VEVs of weak bi-doublets at the GUT-scale contained in 
(E ® S)o.i,27 already discussed in eq. (|37[) follow in a straight-forward manner. 

To have a rough idea of the order of the VEVs involved, using eq. ([37|) and taking the respective Yukawa couplings 
in the range ~ 0.01 - 1.0 and a; ~ 100 GeV, vr = 10" GeV to 10" GeV, Mu ^ 2 x 10^^ GeV, we get (Af) = 10 
MeV to 10 GeV. Our numerical analysis approximately agrees with these results. 



B. Fermion Masses from S'4 Flavor Symmetry 

Investigation on fermion masses and mixings using an S'O(IO) x S'4 model but without any intermediate gauge 
symmetry has been carried out in Ref. [i^ where RG-extrapolated values of charged fermion masses at the GUT scale 
have been used to fix certain model parameters and make predictions in the neutrino sector. Although vr^ — 10^^ — 10^^ 
GeV has been assumed with a view to obtain the right-handed neutrino mass M^ — foVR = 10^^ — 10^"* GeV for 
/o ~ 1, we note that it is difficult to visualise any such value of vr substantially lower than the GUT-scale in a 
single-step breaking scenario; hence the desired value of Mr is not obtainable without adjusting the value of the 
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Majorana coupling to fo 0.001 — 0.01. Also we note that the right choice for the input values of charged fcrmion 
masses is desirable to be at the intermediate seesaw scale rather than the GUT scale. We carry out investigations 
utilising the RG-extrapolated values at Mf{ rj rs 10^^ GeV in the present model where no adjustment of fo is 
needed to obtain the desired sec-saw scale. We find that the experimental data on neutrino mass-squared differences 
and mixings in fact determine the seesaw scale to be Mn ~ 3.78 x 10^"^ GeV. 

Consistent with the 50(10) x ^4 symmetry the superpotential for fermion-Higgs Yukawa interaction is written as, 



W^uk 



+ 4-2*2 + *3*3)(i^)i?0 + /oSo) + ■^(*2*2 " *3*3)(2/l^fl + /l^l) 
+ 4^(-2*i*l + *2*2 + *3*3)(yii?2 + /1S2) 

v6 

+y3[(*2*3 + *3*2)-ff3 + (*1*3 + *3*l)i?4 + (*1*2 + *2*l)-ff5]- 



(39) 



Following the standard notation with Q{Q'") and L{L'~') for left(right)-handed quark and lepton doublets in left-right 

symmetric gauge theory, and denoting Hf and Af as the electroweak bi-doublets in Hi and Sj, we now write the 
Yukawa superpotential consistent with G2213 x S4, just below the GUT-symmetry breaking scale for ^, ~ Mu, before 

the electroweak bi-doublets Af decouple from the superpotential, 



+^T2{yiHt - 3hAt)L^ - LiT2{yiHt - 3/iA^)i^l 
+ -^[-2Qf (yiifl + /iA^)Qf - 2LlT2{y,Ht - 3hAt)L? 

+QlT2{yxHt + hAt)Q^ + L^T2{yiHi - 3/iA^)L^ 
+QlT2{yiHt + hAt)Q^ + LlT2{yM - 3/i At)L^] 

+yz[QlT2HtQ2 + QlT2HtQ^ + QlT2HtQ2 + QjT2HtQ? + QjT2HiQ§ + Q^T2HiQ? 

+{Q ^ L)]. (40) 

The up and down type electro-weak doublets in the six bi-doublets of lO's C 50(10) acquire VEVs < HV^ >= vf 
and < Hf >= vf{i = 0,1, ..,5). The electroweak submultiplets in Ej also acquire induced VEVs < A^ > and 

< Af >,(i = 0,l,2). 

Adding their contributions, the mass matrices of quarks and leptons have the well known forms. 



Ml 



M 



(126) 



3My'^\ = M(io)-3Mii26)^ 



(41) 
(42) 



(43) 



where — fo^R is the degenerate right-handed neutrino mass and we have already noted in Sec. 3 that we can have 
Mjv — Mji substantially below the GUT-scale in this model without having the necessity to adjust the the Majorana 
coupling to be a small fraction of unity . The component mass-matrix elements in the above equations are defined as. 



ao — 2q!2 ^5 (34 

as ao -I- ai -I- a2 as 
a4 a3 ao — ai -|- a2 



(44) 



/3o - 2/32 P5 f34 

Pi P3 P2 



(45) 
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TABLE II: Renormalisation Group extrapolated running masses of quarks and charged leptons of three generations at the 
intermediate scale Mr = 10^'' GeV as estimated in Rcf. 41]. 



tan/3 


10 


55 


TUu (MeV) 


0.888±Ei;l?| 


0.888±o;i79 


rUc (MeV) 


258.094±i:ii 


258.292±i:gl 


mt (GeV) 


94 369±o? il? 

lyTt.uUty -1-25. 833 


104 236±?i ?Si 


TOd (MeV) 


1.829±q;|77 


i.82i±g:|?i 


me (MeV) 


36.426±iiig 


36.289±iin 


mb (GeV) 


1.263±o;og9 


1.576±g:i| 


me (MeV) 


r, OQII 0.0002 

u.^yi±o.ooo2 


n •?SQ-l-0 0005 

u.rf»y±o.ooo2 


m^ (MeV) 


82.553±o;q3| 


82.206±no2 


(GeV) 


1 /ins-i-0.ooo9 

i.4U»±o.oo08 


i.657±8:gl! 



7o - 272 
70 






-71+72 



(46) 



M 



(126) 



Sq - 2S2 

(5o 






61 + 62 

Sa-di+d2 



(47) 



In these equations ai = yi < Hf > , (3i = < Hf >,ji = fi<Af> and 6i = fi < Af > {i not summed). The choice 
of diagonal basis in the down quark sector which automatically also leads to the diagonal basis in the charged lepton 
sector, enables to choose the six parameters, /3i, 5i{i = 0, 1, 2) to be real and = Pi = = 0- AH other parameters 
are, in general, complex. Analytically we express the six real parameters in terms of down-quark and charged lepton 
mass eigen- values at the see-saw scale = Mji), 



f3o 
Pi 
P2 
So 
Si 
S2 



3(mg + m° + m2) + to° + to° + mj^] /12, 



-3m[; + 3m1 - m° + m°] 



3mg + 3m° - Gm^J + m° + ml - 2m°] /24, 
m° + m« + m° - {m° + m« + m^)] /12, 



TOf, -|- m 







2mO 



m" + 2mO] /24. 



(48) 



We utilise the RG-extrapolated values of the running charged fcrmion masses at the intermediate scale fi — Mn « 
vr ~ lO^'^ GeV as shown in Table HIl for tan (3 = 10, 55 [4l[. In the present model the definition tan/3 — v^/vd is valid 
in the presence of MSSM below the intermediate scale. 

Using the down quark and charged lepton masses from Table |TT] and eqs. ([^5]) we obtain. 



449.773 MeV, (3i 



625.971 MeV, P2 = 224.155 MeV, 



(49) 



So = - 15.791 MeV, Si ^ 12.334 MeV, S2 = - 8.074 MeV. 



(50) 



Using low-energy values of CKM matrix elements with its phase S — 60° and using the renormalization factor 
rjv = exp[—{yf^pln{vii/mtop)/l6TT^] — 0.86 leads to the CKM matrix at — Mr — 10^^ GeV, 



Vc 



KM 



0.973852 
-0.227985 - 0.000134610i 
0.00675842 - 0.002851022i 



0.22720 0.00169097- 0.002928801 

0.97301 - 0.000031405^ 0.0369880 
-0.0364054 - 0.000664840? 0.99925 



(51) 
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Defining M„ = diag(m'^, m^, m^) , at first we obtain elements of in terms of the running up-quark masses and 
CKM elements via, 



Mu = V^j^mMuVckm. (52) 

For tan/3 = 10, using eqs.([331), ([?T|l . and ([5^ and Table [Til determines the three parameters ai{i = 3,4,5) 

while three equations are obtained among the other six complex parameters, ai(i = 0, 1, 2) and 7i(i = 0, 1, 2), 

ao + ai+a2+7o + 7i+72 = 369.414±|3;550 _ ,(4.583±i.099)^ 
ao-ai+a2+7o-7i+72 = 9mi M2±lllll-§1 ~ i8. 797 x IQ-"^ , 

ao - 2a2 + 70 - 272 = I7.687±llll ~ t{3m92±''^ill), (53) 



as = -(3423.163±^i:i^)_j(62.685±i:'6?f)^ 
a4 = 634.282±i52^893 _ ■(268.690±6f;225)^ 

as = -{80M8±l'olt'') + ^{9.3M±lill), (54) 
where all parameters are in MeV and the uncertainties in the RHS of these equations reflect the uncertainties in the 



low-energy data 41[. It is clear that the set of three eas. (|53p leaves undetermined three complex (six real ) parameters 
which provide a very rich structure to the model. Because of this, the model may be able to confront the present 
neutrino data and even the future precision data that may emerge from planned and ongoing oscillation experiments. 
On the other hand, it is also possible that the number of parameters may not ensure faithful representation of neutrino 
data because of highly non-linear nature of the problem emerging from see-saw mechanism. 

In order to examine the efficiency of the model in representing the neutrino sector, we use the standard parametriza- 
tion of the Icptonic Pontecorvo-Maki-Nakagawa-Sakata(PMNS) mixing matrix. 



PMNS 



C12C13 S12C13 siae""' 

-C23S12 - S23Sl3Cl2e"^ C23C12 - S23Sl3Sl2e*'' S23C13 
S23S12 - C23Sl3Cl2e"^ -S23C12 - C23Sl3Sl2e*'' C23C13 



diag(e-''^i/2,e-**'^/2,l), (55) 



where Cij = cos0.y, Sy = sin^y, S is the Dirac phase and ipi,ip2 are Major ana phases of neutrinos. These phases have 
range from to 27r. 

We use experimental data on neutino oscillations within the 3(t limit ^42*1: 

0.29 < tan^ 6*12 < 0.64, 
0.49 < tan^ 6*23 < 2.2, 
sin^ 013 < 0.054, 

5.2 < Aml/lQ-^eV^ < 9.8, (56) 
1.4 < AmlJlO-'eV' < 3.4. 

For numerical analysis we exploit the well defined diagonalisation procedure for complex and symmetric mass 
matrices, 

Mi,U* — diag(mi, TO2, TO3), 
U^'M^MlU = diag(m?,TO^,TO^), (57) 

where U is a unitary diagonalising matrix , the light neutrino mass matrix has been defined in ea. (|43p and 
mi{i — 1, 2, 3) are positive mass eigen values. 

For the sake of simplicity we reduce the parameters of the model by treating the parameters 7^(1 = 0, 1, 2) as real. 
Then eq. (|53p determines six real parameters out of a total nine, including real and imaginary parts of ai{i — 0, 1, 2) 
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TABLE III: Fit to the available neutrino oscillation data and predictions of reactor mixing angle 613, leptonic Dirac phase (5), 
Majorana phases (yi, 932) and the CP violation parameter Jcp in the 50(10) x S4 model with see-saw scale at Mn — 3.78 x 10^'' 



GeV and tan /? = 10 


5 (GeV) 


1.025 


1.100 


1.235 


T} (Lre V ) 






A An 
2.4U0 


C (GeV) 




K on 


Zo. lUU 


mi (eV J 


U.OOooD 


U.00o9D 


n AA0A1 
U.OOoOi 


m2 (eV J 


U.009zU 


u.ooyob 


n Ai oc;o 
U.OizDO 


ms (eV) 


U. 05000 


U. 05000 


A A T r* A 

0.07850 




U.UD40 


U.UD/o 


n noooo 
u.uyyzy 




D X iU 


X iU 


y.D X iu 


'-^rriatm (eV ) 


2.5 X 10"'^ 


2.5 X 10"'^ 


3.1 X 10~^ 


sin 012 


0.5io 


n tit G 


A K 1 1 
U.Oii 


sin 023 


0. < is 


U. I io 




sin 013 


0.055 


0.057 


0.052 


5 (radians) 


3.096 


3.048 


3.100 


01 (radians) 


5.67 


5.46 


5.65 


02 (radians) 


5.59 


5.39 


5.65 


Jcp 


2.66 X 10"^ 


6.49 X 10"* 


2.95 X 10"^ 


< rriee > (eV) 


0.00646 


0.00742 


0.00932 


171/3 (eV) 


0.00462 


0.00516 


0.00600 



and real 7^(1 = 0, 1, 2). This choice of parameters implies that the CP-violation has its origin only in the quark sector 
as reflected in the CKM matrix ^] . 

Thus, in addition to the see-saw scale, we are left with three real parameters to fit the neutrino oscillation data 
on four quantities, Atoq, Am^j^^, tan^ 0i2, and tan^ 623 and make predictions on sin ^13, leptonic Dirac phase (S) and 
Majorana phases ( 1, sum of the three light neutrino masses Sm^, the effective matrix element for neutrinoless 
double beta decay, < >, and the kinematic neutrino mass nip to be measured in beta decay where 

< m,e >= , iUpMNs)'^^, ^0 = (E' 1 PpMNsl'm^y^^. (58) 

Equivalently, the three real unknown parameters are defined as, 

^ = 70 - 272, 
V = 70 + 71 + 72, 

C = 7o - 7i + 72- (59) 

Even in the constrained parametrisation of the model , we find that ^, 77 and ^ are quite efficient in describing the 
present neutrino oscillation data. Some examples of our fit to the data and model predictions are shown in Table IIIII 
We find that the see-saw scale is determined to be Mn = 3.78 x 10^'^ GeV for hierarchial neutrino masses. The first 
and the second columns show that for fixed values of rj and C, the parameter ^ is very effective in controlling the value 
of the solar neutrino mixing angle (^12)- Within the uncertainties shown in eq. (|53p and ea. (j54p . the predicted reactor 
mixing angle occurs in the range Ois ~ 3° — 5° which is within the accessible limit of ongoing and planned experiments 



44| . The sum of the three neutrino masses are found to be well within the cosmological bound [45|. The leptonic 
Dirac phase turns out to be closer to tt with S ~ 2.9 — 3.1 radians and the two Majorana phases are within 5.3 — 5.7 
radians. The predicted values of matrix element for double beta decay and the kinematical mass for beta decay are 
found to be nearly two orders smaller than the current experimental bounds 4^, 47, 48|. Similar conclusion has been 
also obtained for hierarchial neutrinos with S4 flavor symmetry in the non-SUSY standard model [28[. The Jarlskog 
invariant [i^ is found to vary between Jcp — 2.95 x 10~^ and Jcp — 10^^ where the smaller (larger) value depends 
upon how much closer (farther) is the Dirac phase (S) from tt. We observe that the predictions of this model in the 
neutrino sector made at the high see-saw scale is to remain stable under radiative corrections when extrapolated to 
low energies especially since the light neutrino mass eigen values are small [50| . 
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VII. SUMMARY AND CONCLUSION 

In this work we have addressed the question of possible existence of R-parity and Parity conserving left-right 
gauge theory as an intermediate symmetry in supersymmetric 50(10) grand unified theory with manifest one-loop 
unification of the gauge couplings. We found that it is possible to have this intermediate gauge symmetry provided 
both the left-right gauge theory and SO(IO) are extended to contain >5'4 flavor symmetry. The particle spectrum needed 
to implement the gauge coupling unification is found to match into different G2213 x >S'4 representations leading to 
exactly vanishing values of two RG coefficients. The Higgs spectrum is also found to be consistent with the mass 
spectra analysis for the SO{10) x ^4 model with fine-tuning of certain parameters in the Higgs superpotential. At first 
ignoring the light neutrino mass constraint, the left-right symmetry breaking scale is allowed to have a wide range 
of values with Mji = 5 x 10^ GeV to 10^^ GeV and no tuning of the Majorana coupling is needed in this model to 
obtain desired value of the see-saw scale substantially below the GUT scale. 

We have carried out analysis of fermion masses and mixings using RG extrapolated values of the low-energy data at 
the intermediate scale ~ 10"'^^ GeV in SUSY SO{10) x S4 for the first time. Even in the case of constrained model 
paramctrisation where CP-violation originates only from the quark sector through CKM matrix, the model is found 
to fit all values of quark and lepton masses and mixings including very large values of mixings and very small values of 
masses in the neutrino sector. The neutrino oscillation data determines the see-saw scale to be Mjv — 3.8 x 10^'^ GeV 
for hierarchial neutrinos. Apart from predictions on leptonic CP-violating parameter, Dirac and Majorana phases, 
the predicted values of the reactor-neutrino mixing angle, ^13 ~ 3° — 5°, are accessible to ongoing and planned long 
baseline experiments on neutrino oscillations. 

It would be interesting to investigate prospects of this model with all complex parameters in type I see-saw and the 
case of experimentally testable quasi-degenerate neutrino spectrum with type II see-saw or a combination of both the 
type I and type II see-saw models in future works. 
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